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MULTI-TRANSONIC BLACK HOLE ACCRETION DISCS WITH ISOTHERMAL STANDING SHOCKS 
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ABSTRACT 

In this work we would like to address the issue of shock formation in black hole accretion discs. We provide a 
generalized two parameter solution scheme for multi-transonic, accretion and wind around Schwarzschild black 
holes, by mainly concentrating on accretion solutions which may contain steady, standing isothermal shocks. Such 
shocks conserve flow temperature by the dissipation of energy at the shock location. We use the vertically inte- 
grated 1 .5 dimensional model to describe the disc structure, where the equations of motion apply to the equatorial 
plane of the central accretor, assuming the flow to be in hydrostatic equilibrium in the transverse direction. Unlike 
previous works in this field, our calculation is not restricted to any particular kind of post-Newtonian gravitational 
potentials, rather we use all available pseudo-Schwarzschild potentials to formulate and solve the equations gov- 
erning the accretion and wind only in terms of flow temperature T and specific angular momentum A of the flow. 
The accretion flow is assumed to be non-dissipative everywhere, except possibly at the shock location, if any. We 
observe that a significant region of parameter space spanned by {A, T} allows shock formation. Our generalized 
formahsm assures that the shock formation is not just an artifact of a particular type of gravitational potential, 
rather inclusion of all available black-hole potentials demonstrates a substantially extended zone of parameter 
space allowing for the possibility of shock formation. We thus arrive at the conclusion that the standing shocks 
are essential ingredients in rotating, advective accretion flows of isothermal fluid around a non-spinning astro- 
physical black hole. We identify all possible shock solutions which may be present in isothermal disc accretion 
and thoroughly study the dependence of various shock parameters on fundamental dynamical variables governing 
the accretion flow for all possible initial boundary conditions. Types of shocks discussed in this paper may appear 
to be 'bright' because of the huge amount of energy dissipation at shock, and the quick removal of such energy to 
maintain the isothermality may power the strong X-ray flares recently observed to be emerged from our Galactic 
centre. The results are discussed in connection to other astrophysical phenomena of related interest, such as the 
QPO behaviour of galactic black hole candidates. 

Subject headings: accretion, accretion disks — black hole physics — hydrodynamics — shock waves — QPO 



1. INTRODUCTION 

Perturbation of various kinds may produce discontinuities in as- 
trophysical fluid flow. Discontinuities in fluid flow are said to 
take place over one or more surfaces when any dynamical and 
/ or thermodynamic quantity may change discontinuously as 
such surfaces are crossed; the corresponding surfaces are called 
surfaces of discontinuity. Certain boundary conditions are to be 
satisfied across such surfaces and according to those conditions, 
surfaces of discontinuities are classified into various categories; 
the most important are being the shock waves or shocks, which, 
when the flow is adiabatic, obey the following conditions (Lan- 
dau & Lifsitz, 1959, LL hereafter): 



[pu]=0, [/:' + pM^]=0, 



^u^ 



■ + w 



:0 



where p,u,p and w are the fluid density, flow velocity, fluid 
pressure and the enthalpy of the flow respectively, and for any 
two arbitrary flow variables a and (3, 

[a + l3] = {(a + l3)_-(a + l3\} 

[a/3] = {a-f3--a+P+} . 

The subscripts (+) and (-) refer to the quantities measured 
just after and just before the flow encounters the surface(s) of 



discontinuity. Quite often, such shock waves are generated 
in various kinds of supersonic astrophysical flows having in- 
trinsic angular momentum, resulting in a flow which becomes 
subsonic. This is because of the fact that the repulsive centrifu- 
gal potential barrier experienced by such flows are sufficiently 
strong to brake the motion and a stationary solution could be 
introduced only through a shock. Rotating, transonic astro- 
physical fluid flows are thus believed to be 'prone' to the shock 
formation phenomena. It has been established in recent years 
that in order to satisfy the inner boundary conditions imposed 
by the event horizon, accretion onto black holes should exhibit 
transonic properties in general; which further indicates that for- 
mation of shock waves are possible in astrophysical fluid flows 
onto galactic and extra-galactic black holes. One also expects 
that shock formation in black hole accretion discs might be a 
general phenomenon because shock waves in rotating astro- 
physical flows are convincingly able to provide an important 
and efficient mechanism for conversion of a significant amount 
of the gravitational energy (available from deep potential wells 
created by these massive compact accretors) into radiation by 
randomizing the directed infall motion of the accreting fluid. 
Hence shocks possibly play an important role in governing the 
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overall dynamical and radiative processes taking place in as- 
trophysical fluids and plasma accreting onto black holes. Thus 
the study of steady, standing, stationary shock waves produced 
in black hole accretion has acquired a very important status in 
recent years. 

The issue of the formation of steady, standing shock waves 
in black hole accretion discs is addressed in two different ways 
in general. First, one can study the formation of Rankine- 
Hugoniot shock waves (RHSW) in a polytropic flow. Radiative 
coohng in this type of shock is quite inefficient. No energy is 
dissipated at the shock and the total specific energy of the ac- 
creting material is a shock-conserved quantity. Entropy is gen- 
erated at the shock and the post-shock flow possesses a higher 
entropy accretion rate than it's pre-shock counterpart. The 
flow changes its temperature permanently at the shock. Higher 
post-shock temperature puffs up the post-shock flow and a 
quasi-spherical, quasi-toroidal centrifugal pressure supported 
boundary layer is formed in the inner region of the accretion 
disc. Fukue (1987) was the pioneer to initiate the systematic 
study of such shock dynamics using a rigorous equation of state 
valid for from the non-relativistic regime to the highly relativis- 
tic regime. Then Chakrabarti and collaborators (Chakrabarti 

1989, 1996 and references therein, Abramowicz & Chakrabarti 

1990, Chakrabarti and Molteni 1993) used the verticaUy inte- 
grated flow model of Matsumoto et. al. (1984) to extended 
Fukue's work towards the global study of shock formation in 
post-Newtonian geometry. However, such shock solutions were 
obtained either on a case by case basis, or, even when successful 
attempts were made to provide a more complete analysis, the 
boundary of the parameter space responsible for shock forma- 
tion was obtained only for global variation of the total specific 
energy £ (or accretion rate M) and specific angular momentum 
A of the flow, and not for variations of the polytropic constant 7 
of the flow; rather, accretion was usually considered to be ultra- 
relativistic', which may not always be a realistic assumption. 
As 7 is expected to have great influence on the radiative prop- 
erties of the flow in general, ignoring the explicit dependence 
of shock solutions on 7 limits claims of generahty. Also, all 
such so-called global shock solutions have been discussed only 
in the context of one parlicular type of pseudo-Schwarzschild 
black hole potential, namely, the Paczyhski & Wiita (1980) po- 
tential ($1 hereafter). For accretion in Kerr geometry, various 
authors (Lu 1985,1986, Lu et al 1987, Nakayama and Fukue 
1989) made attempts to provide a complete general relativistic 
description of shocked accretion flow by making a number of 
assumptions, some of which, however, may not appear to be 
fuUy convincing (see a detail discussion in Das 2002) and it is 
fair to say that no well accepted complete general relativistic 
global shock solution exclusively obtained for a hydrodynamic 
accretion disc around a Schwarzschild BH has yet appeared 
in the hterature. Motivated by these limitations. Das (2002, 
D02 hereafter) provided a generalized formalism allowing one 
to construct and solve the equations governing the hydrody- 
namic black hole accretion of a polytropic fluid in all available 
pseudo-Schwarzschild potentials which may contain steady, 
standing RHSW. In this work, the explicit dependence of shock 
solutions on all important accretion parameters was thoroughly 
studied. The generalized formalism developed by D02 has also 
been used to explain related interesting phenomena hke the 



analytical computation of QPO frequencies (Das 2003) and the 
generation of accretion powered-galactic jets (Das, Rao and 
Vadawale 2003). 

In another class of shock study, one concentrates on shock 
formation in isothermal black hole accretion discs. The char- 
acteristic features of such shocks are quite different from the 
non-dissipative shocks discussed above. In isothermal shocks, 
the accretion flow is allowed to dissipate a part of its energy 
and entropy at the shock surface to keep the post-shock temper- 
ature equal to its pre-shock value. This maintains the vertical 
thickness of the flow exactly the same just before and just af- 
ter the shock is formed. Simultaneous jumps in energy and 
entropy join the pre-shock supersonic flow to its post-shock 
subsonic counterpart. Chakrabarti (1989a, C89 hereafter) stud- 
ied the properties of such shocks in the rotating isothermal 
accretion around compact objects. However, his work was re- 
stricted to only one post-Newtonian black hole potential, i.e., 
the Paczyhski and Wiita (1980) potential ($1) and the depen- 
dence of shock parameters on initial boundary conditions was 
not thoroughly studied. Also, the degeneracy of the shock loca- 
tions was not removed completely. Later on, Yang and Kafatos 
(1995, YK hereafter) investigated the formation of isothermal 
shocks and studied the stability properties of such shocks for 
general relativistic accretion onto Schwarzschild black holes. 
They also, quite successfully, computed the jump conditions 
for relativistic shocks. YK analyzed the stability of shocks in 
a more efficient way compared to the method used by C89 and 
could remove the two-fold degeneracy of shock locations; in 
this way they could properly identify the most stable shock 
location. However, unlike C89, YK used a one-dimensional 
conical model for accretion discs, where the existence of ver- 
tical hydrostatic equilibrium was not taken into account. Also, 
the expUcit dependence of various shock-related quantities on 
accretion parameters was not fuUy explored. Nevertheless, YK 
still deserves special attention because it is the only work in the 
literature which attempts to study isothermal shock formation 
in a fully relativistic frame work. 

In this paper, we would like to address the issue of shock 
formation in isothermal black hole accretion discs from a 
more general point of view. Instead of confining our calcula- 
tion to a particular black hole potential, all available pseudo- 
Schwarzschild potentials will be exhaustively used to investi- 
gate the possibility of shock formation in hydrodynamic, in- 
viscid, isothermal black hole accretion discs. Following Mat- 
sumoto et. al. (1984), we will be using the vertically inte- 
grated 1.5-dimensional model to describe the structure of the 
accretion discs, which is, perhaps, more realistic compared to 
simplified one dimensional conical flow model used by YK. In 
this 1.5 dimensional model (which has extensive use in liter- 
ature, see Chakrabarti 1996 and references therein) the equa- 
tions of motion are written on the equatorial plane of the cen- 
tral accretor whereas the flow remains in hydrostatic equihb- 
rium in the transverse direction. For all black hole potentials 
we first formulate and solve the equations governing multi- 
transonic accretion and wind around the central accretor in 
terms of only two parameters, namely, specific flow angular 
momentum A and flow temperature T. Then we provide a 
self-consistent scheme to capture the stable shock location(s) 
for multi-transonic accretion flow and finally we thoroughly 



' By the term 'ultra-relativistic' and 'purely non-relativistic' we mean a flow with 7=5 and 7=1 respectively, according to the terminology used in Frank et. al. 
1992. 
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study the dependence of various shock related quantities on 
fundamental accretion parameters. In this way we provide a 
generalized two parameter solution scheme for multi-transonic 
isothermal black hole accretion disc for all available pseudo- 
Scharzschild black hole potentials which contains dissipative 
shock waves. Our generahzed formahsm assures that our model 
is not just an artifact of a particular type of potential only, and 
inclusion of every black hole potential allows a substantially 
extended zone of parameter space allowing for the possibility 
of shock formation. 

Rigorous investigation of the complete general relativistic 
multi-transonic black hole accretion disc and wind is believed 
to be extremely complicated. At the same time it is understood 
that as relativistic effects play an important role in the regions 
close to the accreting black hole (where most of the gravita- 
tional potential energy is released), purely Newtonian gravita- 
tional potential (in the form = where Mbh is the 
mass of the compact accretor) cannot be a realistic choice to 
describe transonic black hole accretion in general. To compro- 
mise between the ease of handling of a Newtonian description 
of gravity and the reaUstic situations described by complicated 
general relativistic calculations, a series of 'modified' Newto- 
nian potentials have been introduced to describe the general rel- 
ativistic effects that are most important for accretion disc struc- 
ture around Schwarzschild black holes. Introduction of such 
potentials allows one to investigate the complicated physical 
processes taking place in disc accretion in a semi-Newtonian 
framework by avoiding pure general relativistic calculations so 
that most of the features of spacetime around a compact ob- 
ject are retained and some crucial properties of the analogous 
relativistic solutions of disc structure could be reproduced with 
high accuracy. Hence, those potentials might be designated as 
'pseudo- Schwarzschild' potentials. Our calculations in this pa- 
per will be based on four such pseudo Schwarzschild potentials. 
Explicit forms of those four potentials are the following (see 
Das & Sarkar 2001, D02 and Artemova, Bjornsson & Novikov 
1996, ABN hereafter, and references therein, for detail discus- 
sion about various properties of these potentials): 

2" 



$l(r) = - 



1 



2(r-l) 



$2(r) = 



1 

~2r 



1 



2r 



■ + 12 



$3('') = -l+ 1 



; *4(r)=i/nn-i 



(1) 



where r is the radial co-ordinate scaled in units of 
Schwarzschild radius. Hereafter, we will define the 
Schwarzschild radius rg as 

2GMbh 



(where Mbh is the mass of the black hole, G is universal grav- 
itational constant and c is velocity of light in vacuum) so that 
the marginally bound circular orbit rt, and the last stable cir- 
cular orbit rs take the values 2r^ and 3rg respectively for a 
typical Schwarzschild black hole. Also, total mechanical en- 
ergy per unit mass on r, (sometimes called 'efficiency' e) may 
be computed as -0.057 for this case. We will use a simphfied 
geometric unit throughout this paper where radial distance r is 
scaled in units of rg, radial dynamical velocity u and polytropic 
sound speed a of the flow is scaled in units of c (the velocity of 
light in vacuum), mass m is scaled in units of Mbh and all other 
derived quantities would be scaled accordingly. For simpUcity, 
we will use G = c = 1. Among the above potentials, $i(r) was 



introduced by Paczynski and Wiita (1980) which accurately 
reproduces the positions of and r^. Also the Keplerian dis- 
tribution of angular momentum obtained using this potential 
is exactly the same as that obtained in a pure Schwarzschild 
geometry. $2('') was proposed by Nowak and Wagoner (1991) 
to approximate some of the dominant relativistic effects of the 
accreting black hole (slowly rotating or non rotating) via a mod- 
ified Newtonian potential. It has the correct form of r,, and it 
produces the best approximation for the value of the angular ve- 
locity fij (as measured at infinity) at r, and the radial epicyclic 
frequency k (for r > r^). $3(r) and $4(r) were proposed by 
Artemova, Bjornsson & Novikov 1996, ABN hereafter, to pro- 
duce exactly the same value of the free-fall acceleration of a test 
particle at a given value of r as is obtained for a test particle at 
rest with respect to the Schwarzschild reference frame ($3) and 
to produce the value of the free fall acceleration that is equal to 
the value of the covariant component of the three dimensional 
free-fall acceleration vector of a test particle that is at rest in the 
Schwarzschild reference frame ($4) respectively. Hereafter, 
we will denote any /th potential as ^i(r) where {;= 1,2,3,4} 
corresponds to {^i{r), ^lir), ^sir), ^'4('")} respectively. 

At this point, one may ask why we are interested in study- 
ing shock formation (and the overall structure of the isother- 
mal black hole accretion disc in general) for various pseudo 
Schwarzschild potentials when YK have already performed 
shock study using complete general relativistic frame work. 
The answer is mainly two-fold. Firstly, we have used, perhaps, 
a more realistic flow geometry (1.5 dimensional vertically in- 
tegrated disc model) than that of YK and we have provided a 
more detail dependence of various shock related quantities on 
fundamental accretion parameters. Secondly, and more seri- 
ously, even if someone can provide a completely satisfactory 
model for shock formation in full general relativistic isother- 
mal black hole accretion, still the importance of our work may 
not be under estimated as we believe. Rigorous investigation 
of some of the shock related phenomena is extremely difficult 
(if not completely impossible) to study using full general rel- 
ativistic framework. Hence one is expected to always rely on 
these pseudo-potentials because of the ease of handling them. 
Also, there are possibilities that in future someone may come 
up with a pseudo-Schwarzschild potential better than any of the 
potentials used in this paper, which will be the best approxi- 
mation for complete general relativistic investigation of multi- 
transonic shocked flow. In such case, if one already formulates 
a generalized model for multi-transonic shocked isothermal ac- 
cretion disc for any arbitrary black hole potential, exactly what 
we have done in this paper, then that generalized model will be 
able to readily accommodate that new $(r) without having any 
significant change in the fundamental structure of the formula- 
tion and solution scheme of the model and we need not have to 
worry about providing any new scheme exclusively vaUd only 
for that new potential, if any. 



2. GOVERNING EQUATIONS 

Hereafter we all use the notation ["Pil, for a set of values of 
{A, r} for any particular It is to be mentioned here that 
for spherically symmetric isothermal accretion, one can solve 
the whole problem if only T alone is known. A one parameter 
solution scheme for Bondi (1952) type black hole accretion in 
aU $,s has also been formulated quite recently (Sarkar & Das 
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2002). For isothermal equation of state, 

P='— = c]p ila) 
M 

where P, p, T and Cj are the pressure, density, constant temper- 
ature and the constant isothermal acoustic velocity of the flow 

respectively, R and p being universal gas constant and the mean 
molecular weight, Q and T can be related as 



(2b) 



where 6 = is a constant, niu'^rnp being the mass of the 

hydrogen atom and k is Boltzman's constant. Following stan- 
dard literature, we consider a thin, rotating, axisymmetric, in- 
viscid steady flow in hydrostatic equilibrium in transverse di- 
rection. The assumption of hydrostatic equilibrium is justified 
for a thin flow because for such flows, the infall time scale 
is expected to exceed the local sound crossing time scale in 
the direction transverse to the flow. The flow is also assumed 
to posses considerably large radial velocity which makes the 
flow 'advective' . We also assume that the disc to be non-self- 
gravitating, i.e. the increase of Mbh due to the accretion process 
itself is negligible. The local half-thickness /z,(r) of the disc for 
any ^i{r) can be obtained by balancing the gravitational force 
by pressure gradient and can be expressed as 



I rT 
<I>Kr) 



(3) 



where $,'(r) = ^^jf^- The radial momentum conservation equa- 
tion and the equation of continuity could be written as: 



du dP ^ d 
dr p dr dr 



= 



— [uprhiir)] : 
dr 



m 



where ^"/^{r) is the effective potential which can be defined as 
the sunomation of the gravitational potential and the centrifugal 
potential for matter accreting under the influence of and can 
be expressed as: 



^ff(r) = ^i{r)+--^ (5a) 

However, for inviscid flow, A(r) has been assumed to be A only. 
One can calculate the expression for the effective potential in 
full general relativity (D02): 



■A2(l + r) 



•1, 



(5b) 



to show that unUke ^^/^(r), ^QR(r) cannot be obtained by lin- 
early combining the gravitational and rotational energy contri- 
butions because various energies in general relativity are com- 
bined together to produce non-linearly coupled new terms. Us- 
ing eq. (2a - 2b) and (3), we integrate eq.(4a-4b) to obtain two 
conservation equations for our flow as: 



^ + eTlnpe(r)+^ + $,(r) = C 



Min = Ope(r)Ue(r)r^ 1 



^(r) 



Where C is a constant and M,„ is the mass accretion rate. The 
subscript 'e' indicates the values measured on the equatorial 



plane of the disc; however, we will drop 'e' hereafter if no con- 
fusion arises in doing so and will denote Mi„ by M. 

For a particular value of \V2\v it quite straight forward to 
derive the space gradient of the dynamical flow velocity for any 
as: 



du 
di^ 



36^7- 1 A- 
2r r3 







(7) 



where $" = 4rT^- Since the flow is assumed to be smooth ev- 

' dr- 

erywhere, if the denominator of eq. (7) vanishes at any radial 
distance r, the numerator must also vanish there to maintain the 

continuity of the flow. One therefore arrives at the so called 
'sonic point (alternately, the 'critical point') conditions' by si- 
multaneously making the numerator and denominator of eq. (7) 
equal zero. The sonic point conditions can be expressed as: 



i/=QT-2 = 



A2 



(8) 



where the subscript 's' indicates the quantities are to be mea- 
sured at the sonic point(s). For a fixed ['Pi],, one can solve the 
following polynomial of to obtain the sonic point(s) of the 
flow: 



2A2 



= 



(9) 



Solution of above equation is analogous to solving the follow- 
ing differential equation at critical points x = x^: 



ocy 
where : 



■[{/3Vi(x) + 7V'2(x)} X=Xc +y\ x=xc J y 



f, . . dy(x) .. d^y(x) 



= (10) 



dx ' ■' dx^ ' 

f(x),tpi(x) and ip2(x) are different functions of x and {a,/?, 7} 
are either constant or various functions of initial parameters. 
The critical value of the dynamical velocity gradient at sonic 
point(s) can be computed by solving the foUowing equation for 

(du) . 



1 



du\ 



2 



t 



t 



se^r 3a2 

2r? rj 



(11) 



3. MULTI- TRANSONIC ACCRETION AND WIND 



For all we find a significant region of parameter space 
spanned by [7^2], which allows the multiplicity of sonic points 
for accretion as well as for wind where two real physical inner 
and outer (with respect to the black hole location) X type sonic 
points r,„ and rout encompass one 'O' type unphysical middle 
sonic point r„„v( in between. For a fixed set of {A, T}, one can 
solve eq. (9) to obtain required sonic points and for a particu- 
lar set of values of {A,r}, one obtains at most four real roots, 
one of which is always less than one rg, hence faUs inside the 
event horizon. Thus maximum number of real physical roots 
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is three, hence maximum number of sonic points for any [P2], 
is also three. Apparently there remains an ambiguity in deter- 
mining the values of ['Pi],- for which three sonic points form in 
accretion and the values of [7^2 ], for which three sonic points 
form in wind because just by studying the values of {X,T} or 
{rinjrmidjfout}. One cannot determine whether it is accretion or 
wind which is multi-transonic for that ['Pilr However, we can 
unambiguously classify the parameter space in the following 
way: 

Let us suppose that for a particular value of 17^2 1,, accretion can 
be multi-transonic. So the subsonic far field solution will be- 
come supersonic after crossing the outer sonic point in ac- 
cretion. It would manifest the tendency for passing through the 
inner sonic point rf^'^ again after a shock transition if the solu- 
tion passing through the inner sonic point is more preferred by 
nature for some physical reason(s). One can show that it is the 
total energy content of the flow which will discriminate the so- 
lutions passing through and rf^'^. If the energy content of the 
solution passing through r?^^ is lower than the solution passing 
through r^^, then accretion flow will prefer to pass through r^f^ 
because nature will always prefer a lower energy solution. The 
whole situation would be reversed for multi-transonic winds, 
energy content of wind solutions through rj""'' will be higher 
than that through r^^""'. One can calculate the energy difference 
between solutions passing through rgut and r,„ as: 



Se = eo 



A2 



1 



+ln 



(12) 



where rg = rout and r, = r,„. For a particular if Ai ['Pi],- de- 
notes the universal set representing the entire parameter space 
covering all values of ['Pili, and if B, ['Pi]; represents one par- 
ticular subset of Ai [P2], which contains only the particular val- 
ues of [P2]i for which eq. (9) gives three real physical roots 
(four real roots with one unphysical root lying inside the event 
horizon), then Bi [P2], can further be decomposed into two sub- 
sets d and Vi [P2]/ such that 

CAVil^BAVil only for (5e > 

and: 

C, [P2], CI?, [P2]; only for 5e < 

For [PiJjSCj CPi],, we get multi-transonic accretion and for 
['PiljS'D, \V2\i one obtains multi-transonic wind. In Fig. 1, we 
plot: 

(A„7;)eLP2],eC,LP2JiCi3,[p2j,- 

and in Fig. 2. we plot: 

For all $i, while the specific angular momentum is plotted 
along X axis, the flow temperature T\q, scaled in the unit of 
IQ^^K, is plotted along Y axis. For Fig. 1, AiBjCi represents 
the region C, [^2]; for a particular $, and for Fig. 2, AiBjC;Dj 
represents the region 2?, \V2\i for a particular From Fig. 1 
it is evident that no region of parameter space common to all 
$, is found for which [P2], € C, [P2],. However, significant re- 
gion of parameter space is obtained for which [P2]; G C,- [P2]/ 
for $2 and $3. It is interesting to note that (see Fig. 2 of D02) 
same trend is observed for multi-transonic polytropic accretion 
as well. The maximum flow temperature in 7io for which one 



can obtain multi-transonic accretion may be denoted by Tf^ 
for any ^iir). One observes that 

and the minimum value of angular momentum Af " for which 
one obtains the multi-transonic accretion, can be arranged as: 

Hence one can observe that for $3 and $2. even very low an- 
gular momentum flow (Quasi-Bondi accretion, so to say) may 
produce multi-transonic accretion. 

Situation is not exactly similar for multi-transonic wind. From 
Fig. 2 one observes that there is a common region for 2?, [^2], 
for flows in <l>i, $4 and $3, and there exist a separate common 
zone for 2?, [P2]i in between $2, ^4 and ^j,. Here also same 
sequential order is maintained for T^"'"^: 



7 max ^ rp, 
3 > 



max 
4 



^ rpmax ^ rpmax 



and for Af ": 



ymin ^ ymin ^ ymin 



< XT 



However, there are significant difference for range of A for 
multi-transonic accretion and wind. Unhke multi-transonic ac- 
cretion, three sonic points in wind form for considerably high 
values of A, ; particularly for $1 and ^4, the value of A is very 
large. It is to be noted that for polytropic accretion and wind, 
the situation was quite different (see Fig. 2 of D02); unlike 
isothermal accretion and wind, it was observed there that 



A" 



wind 



accretion 

It is evident from Fig. 2 that multi-transonic wind is usually ob- 
tained obtained only for strongly rotating flows, i.e., for flows 
with intense rotational energy content, and quasi-Bondi flow 
does not produce multi-transonic wind, this trend is more ex- 
plicit for flows in $1 and $4. 

4. SHOCK FORMATION AND RELATED PHENOMENA 

In this work we are not interested in shock formation in 
wind, rather we will concentrate on shock formation in accre- 
tion only. For any if shock forms, then [VtIi responsible 
for shock formation must be somewhere from the region AiBjCi 
(see Fig. 1), i.e, for which [P2]; G C,- [P2];, though not all such 
[P2], will allow shock transition. A multi-transonic accretion 
solution containing shock waves may behave as follows: 

The far-field subsonic solution will become supersonic af- 
ter passing the outer sonic point rout- As the specific energy 
content of the solution passing through the inner sonic point 
r,„ is less compared to the already supersonic solution through 
roui, the solution through r„,„ will tend to 'jump down' to it's 
counter solution passing through r,„. This jump will take place 
via a shock transition and certain amount of energy will be dis- 
sipated at shock, while the flow temperature (hence the accous- 
tic flow velocity) will remain the same. The post shock flow 
will become subsonic and will again accelerate to supersonic 
as it passes through r,„. The combined shocked flow field is 
thus subsonic at infinity and supersonic near the black hole, a 
physicaUy acceptable solution for black hole accretion. 
Like temperature, the baryon flux, the momentum flux and a 
specific function of Mach number are also conserved at shock. 
Following C89, those conservation equation could be written 
as: 

= 0, [P+pu^] =0, 



= 



(13) 
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where M denotes the Mach number, and p and P are the verti- 
cally integrated values of density and pressure. Also 

[a + /3] = {(a + /?)_-(a + /3)J 

[af3] = {a-/3--a+p+} 

as described in §1. Using eq. (13), we can formulate a general- 
ized shock condition, which is the following 

1 u^ + e^T 

P- u_p- + {2p+- p-)<S^T 

The above equation is vaUd only at the shock location, "+" and 
"-" refers to the post and pre shock quantities respectively. For 
the kind of shock discussed here, the shock strength 5, (which 
is basically the ratio of the pre- to post-shock flow Mach num- 
ber) is exactly equal to the shock compression ratio Rcomp (the 
ratio of vertically integrated post- to pre-shock flow density). 



Si=R 



comp 



(15) 

The amount of energy dissipation at shock can be computed as 



Ae = 



{Rlomp-'^)-Q^T In {Rcomp) (16) 



where M is the accretion rate scaled in the unit of Eddington 
rate. In another way eq. (16) can be written in M independent 
form as: 



Ae=i(M2-M^)-e2r/n5,- 



(17) 



Fig. 3 demonstrates few typical flow topologies of the inte- 
gral curves of motion for shocked flow in various (indicated 
in the figure). While the distance from the event horizon of 
the central black hole (scaled in the units of rg and plotted in 
logarithmic unit) is plotted along the X axis, the local Mach 
number of the flow is plotted along the Y axis. One can easily 
obtain such a set of figures for any [Vi\i which allow shock 
transition. For all figures, ABC1C2D represents the transonic 
accretion passing through the outer sonic point rout (marked as 
B) if a shock would not form. However, as 5e in between rout 
and r,„ is greater than zero, the flow will encounter a shock, be- 
comes subsonic and jumps on the branch J2F, which ultimately 
hits the event horizon supersonically after it passes through 
the inner sonic point r„, marked by I in the figure. The small 
circle with a dot at its centre ("0") indicates the location of the 
unphysical middle sonic point r„id. It is to be mentioned that 
by solving eq. (14) one obtains two formal shock locations Ci 
and C2 (vertical lines marked by arrowheads through Ci and 
C2 represent the shock transitions). For every and for any 
\V2\i allowing shock transition, r„id always lies in between two 
shocks passing through C\ and C2. In other words, the 'outer' 
shock at C\ is obtained in between r„„, and rm,j while the 'inner' 
shock at C2 is obtained in rmtd and r,>,. For all <I>,s, it is observed 
that the location of the inner and the outer sonic points anti- 
correlate with A as well as with T. Hotter and strongly rotating 
flows push the real physical sonic points towards the event hori- 
zon. The corresponding values of {r,„, r„„v/, ro„,}, the shock 
locations at C\ and C2 and the corresponding shock strengths 
are indicated at the top of each figure, while the corresponding 
values of A and T (scaled in the unit T\o T /Tiq(= lO^^A")) 
for which the solutions are obtained, are indicated inside each 
figure. GBH represents the 'self wind' of the flow, which, in 
the course of it's motion away from the black hole to infinity, 
becomes supersonic after passing through rout at B. The overaU 



scheme for obtaining the above mentioned integral curves is as 
follows: 

For a particular [V2]i, we first compute r,„, r„„v/ and rout 
by solving eq. (9). Then we obtain the dynamical velocity 
gradient of the flow at sonic points by solving eq. (11). We 
then calculate the local dynamical flow velocity u{r), the lo- 
cal radial Mach number M(r), the fluid density p(r) and any 
other related dynamical or thermodynamic quantities by solv- 
ing the eq. (6-11) from the outer sonic point using fourth 
order Runge-Kutta method. We start integrating from rout in 
two different directions. Along BH, we only solve for uir) and 
M(r), because shock does not form in subsonic flows. Inte- 
gration along BC1C2D involves a different procedure. Along 
BC1C2D, along with our calculation of u{r) and M(r), we also, 
at every integration step (with smallest possible step size), keep 
checking whether eq. (14) is being satisfied at that point. Along 
the curve XJ we calculate p+ and u+. At the same time, along 
the curve BC1C2D, from the initial point whose r coordinate 
is equal to the r coordinate of J, to the final point whose r 
coordinate is equal to r,„, we calculate p— Then we substitute 
these values of p_, p+ and u- in eq. (14) and check for which 
r coordinate the values of p-,p+ and m_ satisfy the equation. 
We then pick up that particular point as our shock location rj/, 
and calculate the value of <S, and Rcomp for the shock location(s). 

Like C89 and YK, we also obtain two real physical formal 
shock locations C\ and Ci- C89 obtained exactly same shock 
strength for all formal shock locations. It is important to note 
that their conclusion is not true. For all and for all values 
of [7^21, allowing shock formation, we obtain distinctively dif- 
ferent shock strength for shocks through Ci and C2. We did not 
follow the stabihty analysis provided by C89 because that kind 
of stabihty analysis cannot remove the degeneracy of shock lo- 
cations, it can only partially remove the multiplicity. We fol- 
low a more robust stability analysis provided by YK where a 
unique stable shock location was obtained. Following YK, we 
find that only the shock through Ci, i.e. the 'outer' shock which 
forms in between rout and r^id is stable, while the other one 
(shock through C2) is not. It is important to note that Nakayama 
(1992) also made important contribution to remove the degen- 
eracy of multiple shocks in relativistic flow. Hereafter, we will 
deal with only the stable 'outer' shock location. In Fig. 4 we 
represents the entire parameter space spanned by [^2], which 
allows shock formation for all The specific flow angular 
momentum A is plotted along the X axis while the flow tem- 
perature T (scaled in the unit of Tio^lO'^/T) is plotted along 
the Y axis. For any pseudo-potential A,B,C, represents the 
parameter space for which [7'2l, G C, {Vj],- The region bounded 
by DiBiCi and marked by S, represents the parameter space for 
which the stable standing shock forms for that particular $, and 
the region A,B,Di represents the region for which accretion is 
multi-transonic yet no stable standing shock forms in accretion. 
If [Pili&Sii'PiliQCiYPili represents the parameter space for 
stable standing shock formation, obviously: 

DiBiCi= [V2I e£i [VilCCi [V2\i (18) 
One can also define [7^2], G.^; [^2],- which is complement of 
£i [Vih related to C; [Vi]; so that: 



AiB,Di 



(19) 



The shock location becomes imaginary in J^i \V2\1 hence no sta- 
ble shock forms in the region A,B,D, even if for [^2], eAjBjD, 
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produces multi-transonic accretion. Hence one should remem- 
ber that although the multiplicity of sonic points in accretion is 
a necessary condition for the existence of shocks, but is not a 
sufficient one. 

We calculate rsi,,Ae,Si,Rcomp and all possible shock related 
quantities for [PiJt GD,B,C, for all <&,s. Hence our calcu- 
lation is much more general one compared to C89 or YK. 
In Fig. 5, we show the variation of r^/, (plotted along the 
Z axis) with A (plotted along the X axis) and T (scaled in 
the unit of Tio and plotted along the Y axis) for all four 

($i^(a), $2->(^), *3^(c), *4-^(^;)). 
One observes that the maximum possible value of shock loca- 
tion is obtained for flows in $i while $3 allows minimum span 
in the value of r^/,. In Fig. (6), we plot the variation of the en- 
ergy dissipation at shock, i.e, the value of Ae (plotted along Z 
axis) with the flow angular momentum A (plotted along the X 
axis) and the flow temperature T (scaled in the unit of Tio and 
plotted along the Y axis) for all four pseudo potentials marked 
in the figure as $1 — >(A),$2 — >(fi),$3 — >(C),<i>4~^(D). The 
shock location correlates with A as well as with T. This is ob- 
vious because a rapidly rotating flow will generate a centrifugal 
barrier at a relatively large distance from the event horizon and 
a hotter flow will generate sufficient amount of thermal pres- 
sure to break the flow well away from the black hole. Also it 
is observed that the stronger the shock (shock with large com- 
pression ratio), the higher the amount of dissipated energy at the 
shock location and that a shock formed closer to the event hori- 
zon will result in a larger amount of energy dissipation. This is 
also obvious because the closer the shock location to the event 
horizon, the higher the amount of available gravitational poten- 
tial energy to be dissipated at the shock. 

5. CONCLUDING REMARKS 

The existing Uterature on shock formation in black hole accre- 
tion discs are heavily inclined towards the study of poly tropic 
accretion and non-dissipative shocks while the formation of 
dissipative shocks in an isothermal flow is a relatively less well 
studied subject. Such a 'preferential investigation' is perhaps 
due to the fact that global isothermaUty in 'black hole accre- 
tion' is difficult to achieve for realistic flows (it is difficult to 
maintain exactly the same temperature throughout a disc), the 
more realistic kind of accretion is expected to be described us- 
ing a poly tropic equation of state. However, one should never 
forget that the techniques developed to study disc geometry and 
shock dynamics in all existing literature (including the present 
paper) are heavily based on various assumptions, and hence, 
are far from being sufficiently complete and self-consistent to 
model the exact physical situation. Even if the vertically in- 
tegrated 1 .5 dimensional disc model developed by Matsumoto 
et. al. (1984), which has been used in this paper, is believed 
to be one of the most successful disc models, it suffers from a 
number of limitations in studying the polytropic shock dynam- 
ics. One of such major limitations is, even for an extremely 
thin shock, i.e. a polytropic shock with effective zero thick- 
ness, the thickness of the flow (i.e. the disc thickness) changes 
abruptly at the shock (pre and post shock half thickness of the 
polytropic flow is quite different) and it is impossible to write 
down the exact pressure balance condition at the shock by ver- 
tical averaging. On the other hand, for isothermal flows with 
dissipative shocks (the kind of shock dynamics studied in this 
paper) at constant acoustic velocity, the flow thickness remains 
the same (because of the quick removal of dissipated energy) 



and a correct pressure balance condition is obtained. Hence we 
believe that the study of shock formation in isothermal flow is 
quite important in its own right and we have provided /or the 
first time a generalized formalism to study such shocks for all 
available pseudo-Schwarzschild black hole potentials. Also we 
have studied the complete disc dynamics for multi-transonic 
accretion and wind and have demonstrated all the shock so- 
lutions and successfully studied the dependence of important 
shock quantities on every fundamental accretion parameter. 
However, in this work is has not been possible to compare the 
different potentials in connection to the suitability of mimick- 
ing the exact general relativistic solutions as was done in D02. 
It was beyond the scope of this paper to formulate and solve 
the complete general relativistic set of equations governing the 
isothermal accretion and wind, we hope to present that work 
somewhere else. 

It is to be mentioned here that a number of recent works 
study the dynamical and radiative properties of the post-shock 
region of advective accretion discs (inner part of the disc, so 
to say) in order to understand the spectral properties of the 
black hole candidate (Shrader & Titarchuk 1998, and refer- 
ences therein) and to theoretically explain a number of diverse 
phenomena which are supposed to take place in the inner part 
of the shocked black hole accretion discs. (Titarchuk, Lapidus 
& Muslimov 1998, and references therein. Das 1998, Das & 
Chakrabarti, 1999, Das 2003, Das, Rao & Vadawale, 2003, 
Naik et.al. 2001). The importance of studying the physics 
of post-shock accretion flow has also been realized from ob- 
servational evidence (Rutledge et. al. 1999, Muno, Morgan 
& Remillard 1999, Webb & Malkan 2000, Smith, Heindl and 
Swank 2002). Thus we believe that our generalized model 
presented in this paper may have far reaching consequences in 
handhng various shock related phenomena because this model 
deals with all available black hole potentials. Nevertheless, one 
limitation of our model is that we have assumed the flow to 
be totally inviscid and viscous dissipation of azimuthal angu- 
lar momentum was not taken into account. Exact modeling of 
viscous multi-transonic black hole accretion, including proper 
heating and cooling mechanisms, is quite an arduous task, and 
we did not dare to attempt that in this paper. However, qual- 
itative calculations show that the introduction of viscosity via 
a radius-dependent power-law distribution for angular momen- 
tum only weakens the strength of the centrifugal barrier and 
pushes the shock location closer to the event horizon, keeping 
the overall basic physics concerning the shock dynamics and 
related issues unaltered; details of this work is in progress and 
will be discussed elsewhere. 

In our work, some of the examples of the shocked multi- 
transonic flow, especially for flows in $2 and $3, contain low 
intrinsic angular momentum. Such weakly rotating flows may 
be found in nature in systems Uke detached or semi-detached 
binaries fed by accretion from OB stellar winds (Illarianov & 
Sunyaev 1975, Liang & Nolan 1984), semi-detached low-mass 
non-magnetic binaries (Bisikalo et al. 1998) or super-massive 
black holes fed by accretion from slowly rotating central steUar 
clusters (Illarionov 1988, Ho 1999 and references therein). 

It is now a well established fact that a number of galactic 
black hole candidate show QPO behaviours. Associated large 
amplitudes and short time scales of such intensity variation 
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suggests that such oscillations are a diagnostic of the accretion 
processes in the innermost regions of the accretion discs of the 
central accretor (van derkUss, 1989, and references therein) and 
it has been suggested in recent years that such QPOs may be 
generated because of the shock oscillation in black hole accre- 
tion discs (Molteni, Sponholz & Chakrabarti 1996, Titarchuk, 
Lapidus & Muslimov 1999) We apply our shock solutions ob- 
tained in this paper to compute the frequency of the QPOs 
generated by such shock oscillations (Das 2003) and study the 
dependence of such frequencies on various fundamental accre- 
tion and shock parameters. 

We would like to complete our discussion by pointing out 
one important finding in our work. Since the accretion flows 
discussed in this paper are chosen to be inviscid and non- 
dissipative, most of the binding energy of the flow is dissi- 
pated at the shocks. Thus such shocks can serve as the efficient 
provider of effective radiative cooling mechanism in the accre- 
tion discs. From Fig. 6 it is evident that for some region of 
parameter space spanned by [7^2]!, the amount of energy dissi- 
pated at the shock could be as high as more than ten percent 
of the total rest mass energy of the flow. This large amount of 
energy has to be immediately removed from the accretion disc 
(at least from the equatorial plane) to maintain the isothermal- 
ity of the flow. This dissipated energy may be drained from the 
disc through various processes Uke a sudden X-ray burst (flare) 
of energy from the accretion discs close to the event horizon 
(few Vg away) or through energetic leptonic/ baryonic outflows 



and jets, the formation of which is likely to be initiated by this 
huge amount of dissipated energy at the shock. Observational 
signature of such flares from the close vicinity of our Galac- 
tic centre black hole has already been reported quite recently 
(Baganoff et al. 2001). The observed length scale of the ac- 
cretion flow from where such flares are generated, are in fairly 
close agreements with our theoretical calculations presented in 
this paper. Thus one concludes that the kind of shock discussed 
in this paper could be efficiently radiative and 'bright'. The 
energy dissipation at the shock and it's channeling from the ac- 
cretion disc system may explain some interesting related issues 
like shock generated X-ray flares and the QPO of galactic mi- 
croquasars (Das & Rao 2003, in preparation) and the genera- 
tion mechanism of accretion powered galactic jet (Das 2003a, 
in preparation). 
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Fig. 1: Parameter space division for multi-transonic accretion in four different pseudo-Schwrzschild potentials marked in the figure. 
Specific angular momentum of the flow (A) is plotted along the X axis while the flow temperature is plotted along the Y axis in the 
unit of lO^'^K. Note that there is no common set of {A, T} for which multi-transonic flow occurs for all potentials, see text for details. 
Fig. 2: Multi-transonic isothermal wind as a function of { A, T — >T / Tio (= IO^'/T) } for various black hole potentials marked on the 
figure. Note how, unUke accretion, multi-transonicity is obtained in wind for very strongly rotating flows (flow with large value of 
A). 

Fig. 3: Solution topology of multi- transonic shocked accretion along with it's 'self- wind' for four different black hole potentials 
marked in the figure. Among the two formal shock solution (marked by vertical lines with downward arrow having different values 
of shock strength, solution through Ci comes out to be the stable one, see text for detail. 

Fig. 4: Global parameter space for accretion solution with shocks (region marked by S\ and bounded by BiCiDi for the potential $i 

and so on) embedded in region for general multi-transonic accretion solutions (region bounded by AiBiCi for <I>i and so on) for four 
different potentials. This figure demonstrates all possible {A, T — >T /T\q (= IO'^/T) } for which shock may form in isothermal black 
hole accretion disc. 

Fig. 5: Variation of shock location (scaled in the unit of Schwarzschild radius and plotted along the Z axis) with specific angular 
momentum (plotted along the X) axis and flow temperature T (scaled in the unit of T\q = lO^^A" and plotted along the Y axis) for four 
potentials $i(a), ^lib), 4>3(c) and ^4,{d)- See text for detail. This figure is complementary to Fig. 4 in the sense that shock locations 
in this figure are obtained for all { A, T — >T / T\o {= IO'^TT) } in Fig. 4 for which shock may form. 

Fig. 6: Variation of the amount of energy dissipation Ae at the shock location (plotted along the Z axis) as a function A (X axis) and 
T/TiQ (= lO'O/:) Y axis for four different potentials $i(A), $2(5), $3(0 and $4(£'), see text for detail. 
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